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ABSTRACT 



. s i 

This thesis lays the foundation for the dynamic stall 
Investigation being conducted at the Fluid Mechanics Lab- 
oratory at NASA-Ames Research Center. Using existing optical 
and electrical equipment, a new dedicated Mlcro-VAX computer 
and Labstar software, an Indraft transonic wind tunnel and 
able technicians to make the proper Interface hardware, the 
project came together In a new test facility at the Fluid 
Mechanics Laboratory. The goal of the thesis was to obtain 
both qualitative and quantitative Information about the wake 
profiles of an airfoil In steady state operations at varying 
angles of attack and tunnel conditions. To accomplish this 
task, schlleren photography was used to obtain a qualitative 
picture of the flow field. With this Information, a two com- 
ponent Laser Doppler Veloclmeter was set up to accurately 
measure the velocity profiles that correspond to the schlle- 
ren photographs. Once this preliminary work Is completed, 
the same apparatus will be used to further Investigate the 
unsteady dynamic stall phenomenon. 
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I. INTRODUCTION 



A. DESCRIPTION OF THE OVERALL DYNAMIC STALL PROJECT 

Static stall Is a well understood and documented phenom- 
enon. A wing or rotor blade that Is operating at a high a, 
angle of attack, will experience a point where separated 
flow will cause an accelerated decrease In lift and a large 
rapid increase In drag. This a at which the wing will pro- 
duce the maximum lift, C] max , Is a very Important aspect 
of airfoil performance. This stall a determines the stall 
speed of the aircraft. Considerable effort has been expended 
to find ways to decrease the stall speed by Increasing the 
Cimax* The majority of this effort has been directed to- 
ward a gradual approach to the stall a, In other words, 
static stall research. 

The wing can be made to continue to produce lift past 
the static stall a If the wing angle of attack Is In- 
creased rapidly past the stall a. The phenomenon of 
dynamic stall Is manifested In the rapid decrease of the 
lift at the higher stall a. There are four basic stages of 
dynamic stall that have been documented [Ref. 13: 

1. Tra 1 1 1 ng-edge flow reversal, progressing forward along 
the a 1 rf ol 1 . 

2. Formation, growth and shedding of the vortex. 

3 . Full sta 1 1 . 

4. Boundary layer reattachment. 

The understanding of this phenomenon will be beneficial to 
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the helicopter, where constantly oscillating rotor blades 
are consistently In the dynamic stall region throughout 
their revolution, and jet fighters, where agile and rapid 
maneuvering through the stall a Is required for surviva- 
bility. By quantifying the velocity field around and behind 
the airfoil, a better understanding of the dynamics of this 
phenomenon will be gained. This new analytical Information 
will provide valuable Information to the design engineer to 
expand the flight envelope of future fighter aircraft and 
improve the performance of the helicopter rotor system. 

B. THESIS GOALS 

This thesis lays the foundation for the dynamic stall 
Investigation project being conducted at the Fluid Mechanics 
Laboratory at NASA-Ames Research Center. Using existing 
optical and electrical equipment, a new dedicated Micro-VAX 
computer, new Labstar software, an existing transonic in- 
draft wind tunnel and able technicians to make the proper 
interfaces, the project came together in a new test cell of 
the Fluid Mechanics Laboratory. The goal of this thesis was 
to first obtain stroboscopic schlieren photographs of a sta- 
tionary airfoil at different a's to obtain qualitative 
data on the flow past the airfoil and the size, shape and 
extent of the wake and then set up a Laser Doppler Veloci- 
meter to obtain quantitative velocity data in order to con- 
struct an accurate velocity profile of the flow field. 



C. SCHLIEREN AND LASER DOPPLER VELOCIMETRY BASICS 

i . SghI leren 

Many optical flow visualization techniques have been 
developed for gas flows. Each of these techniques are sensi- 
tive to variations In density and derivatives of density In 
the flow. Interferometry Is sensitive to the density of the 
gas flow; f. Schlieren Is sensitive to the density grad- 
ient or the change In density with respect to distance; 

y. Shadowgraphy Is sensitive to the second derivative 
of the density; ft] ] the techniques rely on 
the fact that the speed of light Is a function of the local 
Index of refraction In that medium. In gases, the Index of 
refraction Is a function of density. 

The schlieren, German for streak, technique allows 
one to view, qualitatively, the density gradient In a given 
flow field. This visualization system Is considerably sim- 
pler and cheaper to use than other techniques. Therefore, It 
was chosen as the preliminary technique to obtain qualita- 
tive Information about the air flow around and In the wake 
of the airfoil. The schlieren setup Involves using a line 
source of light and projecting It through the test section 
where the flow field alters the path of the light. The line 
source Is usually a monochromatic gas discharge tube. The 
light is projected through a narrow slit to provide a line 
of light of uniform illumination. The dimensions of the slit 
control the definition of the Image seen on the screen. A 
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high quality concave mirror Is usually used to direct the 
light from the source through the test section. On the other 
side of the test section, another similar mirror Is used to 
reflect the light to a viewing screen or camera. Each point 
of the line source will illuminate each point In the test 
section if the source Is collimated and of equal optical 
distance on both sides of the test section. The light from 
the second mirror Is focused onto a knife edge. The place- 
ment and adjustment of the knife edge Is very critical as 
the sensitivity of the entire system depends on It. The 
knife edge Is placed In the plane of the source Image and 
Interferes with a portion of this Image. A portion of the 
light from each point In the test section Is blocked which 
decreases the light Impinging on the viewing screen. The 
density gradients encountered In the test section normal to 
the knife edge cause a refraction of light allowing more 
light from that portion of the test section to pass over the 
knife edge and Illuminate the viewing screen. By rotating 
both the line source and the knife edge both directional 
derivatives of density can be viewed. Also, through the 
lateral adjustment of this knife edge at the source image, 
a fine control over the degree of detail resulting from this 
technique can be exercised. 

The viewing screen can be any Image recording 
device. Polaroid cameras can be set up for Instantaneous 
results or 35 mm cameras can be mounted at the focal length 
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of the lens to obtain numerous photographs during the 
testing. Video cameras can also be used. 

2. Laser Doppler Veloclmetrv 

The techniques for measuring velocity In the past 
used sensors, which need to be Introduced into the flow. 
This necessarily meant Intruding Into and affecting the 
flow. During the last two decades, a new technique, Laser 
Doppler Veloclmetry or LDV, has been available. LDV's are 
Indeed a better technique. It Is nonlntruslve and Is capable 
of directly measuring velocities of particles In the flow. 
Laser Doppler Veloclmetry has been very successful In meas- 
uring flow velocities In turbulent and reverse flows. It has 
been successfully used for the past two decades to obtain 
velocity measurements accurately. 

Laser Doppler Veloclmetry uses two coherent* light 
beams from the same source, which Intersect to form the 
measurement volume. The two light beams form an Interference 
fringe system that Is crossed by particles moving with the 
flow. The light scattered by the particles Is Doppler 
shifted^. The frequency of this scattered light Is 
directly proportional to the particle velocity. With 



1. Coherent Light = light In which the electromagnetic 
waves maintain a fixed phase relationship. (Random House 
College Dictionary) 

2. Doppler Shift = the magnitude of the apparent change 
In frequency resulting from the relative motion of the 
light source and the receiver. (Random House College Diet) 
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suitable selection of flow particles, known also as seeds, 
the fluid velocities can be measured from the particle 
velocities. Using photodetectors, the light scattered by the 
particles is converted into electrical signals. This signal, 
after processing, provides frequency information from which 
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II. EXPERIMENTAL APPARATUS AND PROCEDURE 



A. SCHLIEREN SETUP 

The schlieren system for this experiment was set up 
using both a continuous light source and a stroboscopic 
source. The continuous source was used to provide a time 
averaged schlieren photograph and also to allow use of a 
video tape recorder to record the time history of the flow. 
To provide Instantaneous schlieren photographs of the flow 
field, the stroboscope was used. See Figure 2.1 for the 
schematic of the setup. 

To obtain the 35 mm photographs of the flow field, a US 
Scientific Strobrlte Stroboscope was used. The strobe had 
five settings, 20 Hz to 1 KHz, to adjust the duration of the 
strobe. The 70 Hz setting proved to be the optimal setting. 
A mechanically adjustable slit In front of the strobe was 
used to produce a line source of light needed for the tech- 
nique to operate effectively. A large concave mirror with a 
305 cm focal length reflects this light sheet towards the 
test section. The test section had optical access through 
both sides of the tunnel. The test section's dimensions were 
13 cm by 25 cm and was constructed out of 2.5 cm thick 
Plexiglas. The light passing through the test section was 
then reflected by an Identical mirror located opposite the 
first mirror on the other side of the test section as shown 
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Figure 2.1 Schlleren Schematic 

In Figure 2.1, Light from the second mirror was focussed on 
a vertical knife edge whose position could be adjusted to 
get a highly sensitive system. The knife edge was oriented 
vertically to obtain the visual signature of the density 
gradient in the vertical plane of the test section. After 
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the knife edge, a 200 mm lens was. used to focus the Image 
onto the photographic device. 

Since the test section Is 25 cm wide, the photographs 
obtained from the system actually represent density gradient 
effects averaged over the test section width, as Is always 
the case with the schlleren technique. 

A Nikon 35 mm camera body with motordrive (but no lens) 
was used to record the schlleren Images. The camera was 
placed at the focal point of the 200 mm lens. Fine adjust- 
ments to focus the Image In the camera were made by a micro- 
meter slide on the camera mounting mechanism. Pictures were 
taken by opening the shutter with a remote shutter release 
cable, before the strobe was fired remotely and then closing 
the shutter. With a little practice, this procedure was 
quite rapid. Photographs were taken at two different Mach 
numbers, 0.30 and 0.45, and at three different s; 0°, 
7.7°, 15.5°. See Figures A.l through A. 6, Appendix A, 
for examples of schlleren photographs. This series of photo- 
graphs was studied to analyze the flow qualitatively and to 
determine the locations In the wake for obtaining the veloc- 
ity data for the two different Mach numbers. 

B. LASER DOPPLER VELOCIMETER SETUP 

This section will discuss the specific LDV system setup 
In test cell number four at the Fluid Mechanics Laboratory. 
The following sections will be discussed: transmitting 
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optics, receiving optics and photodetectors, test section, 
signal processing and data processing. Figure 2.2 shows a 
photograph of the transmitting optics and Figure 2.3 shows a 
photograph of the receiving side. Individual schematics will 
be discussed in detail within the following sections. 

The laser used is a Spectra Physics, Model 164, 5 Watt 
Argon-Ion Laser capable of outputting light from 0.351 to 
0.528 microns and up to twelve discrete wavelengths simul- 
taneously CRef. 2 : p . 233. By choosing the blue beam, 0.4880 
microns, and the green beam, 0.5145 microns, for the two 
component Laser Doppler Velocimetry System, nearly 75% of 
the initial output power could be used. This standard 
practice was adapted in this case. 

1 . Transmitting Optics 

Owing to traversing requirements, the laser was 
higher on the optics table than the 10.8 cm optical height 
of the transmitting optics. A set of two mirrors at 45° 
angles to the direction of the incident light were used to 
lower the beam to the require height. Figure 2.4 shows a de- 
tailed schematic of the entire transmitting optics system. 

The vertically polarized beam exiting the multiline 
laser, after being directed by mirrors, is passed through a 
beam collimator, part of the TSI Model 9106 set. Incident 
light is collimated if the beam diameter does not change as 
it passes through space. In reality, however, it is suffic- 
ient if the collimator can Insure this quality through the 
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Figure 2.2 Transmitting Optics Photograph 




Figure 2.3 Receiving Optics Photograph 
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Figure 2.4 Transmitting Optics Schematic 

relatively short optical path of the system. Proper col 11- 
matlon is said to be achieved when the waist* of the beam 
coincides with the focal point of the transmitting lens. By 
sliding the collimator sleeve in and out and observing the 
beam "spots" on a distant wall, the collimator can be set. 



1. Waist = the point on a beam where the diameter Is 
the smal 1 est . 
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From the collimator, the beam passes through a po- 
larization rotator and an attenuator. The polarization ro- 
tator changes the polarization of the Incoming laser beam to 
conform to the polarization requirements of the next optical 
component, a dispersion prism. The attenuator works In 
tandem with the rotator to Insure passage of only the horiz- 
ontally polarized light through these components. The dis- 
persion prism separates the wavelengths, or colors, of the 
light. The two primary colors of this laser, green (0.5145 
microns) and blue (0.4880 microns), are used to form the two 
component system. A set of mirrors, TSI Model 9107, Is used 
to direct each of the beams to the appropriate axis for the 
transmitting optics as shown In Figure 2.4. As can be seen 
from the figure, the green beam Is directed down the center- 
line of the optics and the blue Is directed 25 mm to the 
right of the green. 

The centerline, green, beam passes through a polar- 
ization rotator, TSI Model 9102, to rotate the beam polari- 
zation 90°, then a beam splitter, TSI Model 9115. The 
polarization change Is necessary to Insure proper splitting 
Into two beams from the Incident light. The beam splitter 
divides the Incident beam Into two beams of equal Intensity 
and parallel to the optical axis, spaced 50 mm apart. 
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Since the flow field to be Investigated will have a 
turbulent wake with possible separated flow regions, a freq- 
uency shift system is used for each of the two colored 
beams. Upshifting of the frequency of the laser light prior 
to entering the test section will allow detection of the 
direction the particle traverses through the probe volume. 
The shift results In a moving fringe pattern that will gen- 
erate a photodetector signal equal to the shift frequency 
when measuring a particle with zero velocity. This system 
allows a greater range of measurements of negative veloc- 
ities and, thus, is required for this system to obtain valid 
measurements . 

The pair of green beams now enter a frequency shift- 
ing device called the Bragg cell. This device is an acousto- 
optic unit that upshifts the frequency of an incident beam. 
The shift is set at 40 MHz. This device is placed in the 
optical path of one of the two beams forming the probe vol- 
ume for that color. The outgoing beam will be shifted by a 
40 MHz shift relative to the other beam. The Bragg cell 
operates by deflecting one beam In a fan-like pattern in 
the plane intersecting both incident beams. Each separate 
deflected beam has a frequency that is a multiple of 40 MHz 
from the incident beam. The Bragg cell should be adjusted 
to shift the light intensity to the first order beam. This 
is accomplished by tilting the Bragg cell in its mount. 
This first order beam should have an intensity greater than 
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90% of the Incident light beam [Ref. 3:p. 1483. An optical 
wedge Is Included In this device at the output end of the 
Bragg Cell to steer the deflected, first order beam parallel 
to the unshlfted beam. The result of this Bragg shifting Is 
an Interference pattern In the probe volume that will move 
relative to a fixed particle. This moving fringe pattern can 
be used to differentiate the direction of motion through the 
probe volume. Particles moving with the direction of motion 
of the fringes gives a lower frequency and against the 
fringe movement gives a higher frequency as shown in the 
following equation. The particle moving through the probe 
volume will produce a signal whose frequency Is 

f slgnal = lf shlft +/ " f Doppler‘ 

where the sign is + for particles moving in the opposite 
direction to the fringe movement [Ref. 2:p. 673. 

For optic tube support purposes a TSI , Model 9176, 
Support is used. This device also serves to increase the 
effectiveness of the beam stops of the blue beam by inc- 
reasing the deflected beam separation distance. 

The off-axis blue beam is now brought to the center- 
line using a beam displacer, TSI Model 9174. This model em- 
ploys a rhomboid prism. It displaces the beam a fixed, 25 
mm, distance to the centerline of the optics while keeping 
the beam parallel to the axis. The blue beam then passes 
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through the same series of optical components as the green 
beams: polarization rotator, beam splitter and Bragg cell. 

The next two elements in the optics train are beam 
steering modules, TSI Model 9175, one for each of the 
colors. These modules allow for steering of the shifted beam 
to Insure proper beams crossing. 

The Bragg cell produces several beams, of which 
only the first harmonic, the beam that is shifted by 40 
MHz, is Important. All the other unwanted beams are blocked 
off by means of a beam stop. In the present case, a beam 
stop was used on each of the two shifted beams. 

Following the beam stops, beam spacers, TSI Model 
9113-13, are used to for each set of beams to reduce the 
distance between the beam pairs forming the probe volume. 
The spacing was reduced to 13 mm by using these prisms; 
while maintaining polarity, coherence, beam diameter and 
parallel optical paths. The primary use of these elements is 
to control the fringe spacing in the probe volume by 
changing the angle of intersection of the beam. 
Additionally, however, they permit transmission of the beams 
into the test section by requiring only small mirrors and 
optic components. 

The four beams from the optic tube are directed to 
the focusing lens with five mirrors. The first two move the 
beams from 10.8 cm to 23 cm above the optical table and the 
next two reflect the beams back 180° toward the last 
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mirror in this series. This mirror directs the four beams 
vertically toward the transmitting optics lens. The four 
beams then enter the transmitting optics lens to focus the 
four beams to a point. A final mirror is used to direct this 
focused beam to the test section. The focusing lens and the 
last three mirrors are mounted on a traversing mechanism to 
enable traversing of the focused beams in any of the three 
orthogonal directions as shown in Figure 2.6. 

2. Receiving Optics and Photodetector 

All the receiving optics are mounted on a traversing 
mechanism which is slaved to the traversing mechanism on the 
transmitting side. The two sides were aligned to within 0.01 
mm and could be moved together. A Taskmaster controller is 
used to move the two traverses both manually and through the 
use of computer programs. See Figure 2.5. 

The receiving optics are used to capture the light 
scattered by the particles crossing the probe volume and 
focus this light onto photodetectors. It consists of several 
components; that are mounted on a hardware arm connected to 
the traversing mechanism. The scattered light is collected 
by a series of three lenses that collimate the light and 
direct it toward the polarization color separator. The last 
lens acts as a spatial filter in addition to being a beam 
collimator. Its function is to allow light focused on the 
lens to pass through while blocking all the rest of the 
light. It is used to reduce the stray light scattered from 
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around the probe volume, as Is the case in measurements 



close to a solid surface. 



The polarization beam separator separates the two 



wavelengths present In the scattered light into its con- 



stituent colors; green and blue. This splitter reflects 



light with horizontal polarization upward while transmitting 



the vertically polarized light through It. This insures sep- 



aration of the two primary colors in the scattered light. 



Two notch filters are used at this stage, on each of 



the colors to eliminate residual color contamination. The 
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Figure 2.5 Receiving Optics Schematic 





separated scattered light Is now focused on to the photo- 
detector aperture with a lens. 

The photodetector converts the light signal Into an 
electric signal. The photodetector used Is a photomultiplier 
tube or PMT. The advantages of the PMT over the other types 
are low noise, high gain and broad bandwidth while Its 
fragility and the need for a high voltage source are Its 
disadvantages [Ref. 2:p. 453. 

3 . Test Section 

The transonic test section is formed by Inserting 
upper and lower ventilated walls In a subsonic Indraft wind 
tunnel. The test sectlon / s dimensions are 13 cm by 25 cm. 




Figure 2.6 Transmitting Side of Test Section 
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A 7.5 cm chord NACA 0012 airfoil was mounted between the 
side walls at the centerline. This airfoil Is attached at 
quarter chord to a external device to oscillate the airfoil 
through various angles of attack. The side walls of the test 
section are made of 2.5 cm thick plexlglas to allow complete 
and unobstructed optical access to the airfoil. See Figure 
2.6. The mechanism to mount the model partially blocks the 
views of the leading edge of the model. However, this was 
not a limitation for the present study. 

4. Signal Processing 

The signal from the photodetector must be processed 
to obtain the frequency Information In the scattered light. 
Figure 2.7 shows the block diagram of the signal processing 
components. The signal from the PMT Is amplified by an EIN 
Model 403LA, linear amplifier to boost signal strength by 37 
dB to transmit the signal over a long distance and to be of 
appropriate strength for signal processing. 

A downmixing circuit Is used next to decrease the 
effective frequency shift from 40 MHz to around 10 Mhz to 
Improve the counter resolution. The two elements of the 
downmlxer are a Tektronix S6503 Levelled Sine Wave Generator 
and a Hewlett-Packard, Model 10534A, Passive Diode Mixer. 
The expected Doppler frequency should be estimated to prop- 
erly set the mixer frequency. Using the range of Mach num- 
bers expected as -0.2 to +0.5, the free stream velocity 
range can be calculated using the set of equations In 
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section III A and 20° C as the ambient total temper- 
ature. The velocity range resulting from this calculation Is 
-67 to +167 meters per second. With the fringe spacing of 
19.52 microns and the second equation on page 35, the 
estimated Doppler frequencies are -3.4 to +8.6 MHz. Using 
the equation on page 34, the signal and mixing frequency 
should be expected to fall within the range of 36.6 to 48.6 
MHz. By using a starting mixing frequency of 35 MHz, the 
calculated signal frequency reaching the counters will 

Signal and Data Processing 




Figure 2.7 Signal and Data Processors 
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be 1.6 to 13.6 MHz. The mixing frequency should be adjusted 
whenever the tunnel conditions change or In attempts to 
obtain better counter resolution. 

The next signal processing unit encountered is the 
Series 3000 Macrodyne LDV Processor. This device has several 
functions: variable high and low pass filters. Input and 
output ports, threshold and overload adjustments and freq- 
uency counter with fringe passing selection switch. Both the 
low pass filters, to remove the high frequency noise, and 
the high pass filters, to reduce the low frequency pedestal 
from the signal, are adjustable to provide the optimal fil- 
tering of the signal. Once the signal has been thus condi- 
tioned, a frequency counter is used to measure the signal 
frequency and convert it to digital time data that can be 
sent to the data processor. The counter operates on two 
different systems to validate the Doppler signal: a three 
level detection scheme and a 5/8 or 10/16 comparison. The 
three level detection scheme shown in Figure 2.8 must be met 
before a signal is validated. A threshold level is set on 
the counter to be Just above the level of the steady state 
noise. The signal must pass the following gates to be ac- 
cepted: (1) exceeds the threshold setting <2> crosses zero 
and (3) exceeds a negative threshold. To pass this check, 
the selected number of fringe crossings must sequence 
through each check without error. This detection scheme is 
notably effective in noise discrimination. Once a signal is 
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received and passes the first check, the counter starts to 
measure the Doppler period by opening a time gate and stop- 
ping at the switch set number of cycles. This setting of 
cycles Is able to be selected from a switch on the counter 
control panel. Validated data from the switch being In the 
10/16 cycles position tend to be better for strong Doppler 
bursts because the greatest number of cycles measured tends 
to reduce the effects of noise in the signal. However, this 
setting also restricts the motion of the particle through 
the probe volume and thus Introduces artificial biases of 
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Figure 2.8 Three Level Counter Detection Scheme 
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the data. The particle must traverse the probe volume at 
Its largest length to Insure validated data for 16 fringe 
crossings. For this reason, present experiments were carried 
out with the switch set to the 5/8 position. 

5. Data Processing Hardware 

Once the data Is received from the Macrodyne signal 
processors, It must be multiplexed and transferred to the 
computer for velocity calculations. The NASA LDV Computer 
Interface Is now used to control the multiplexing of data, 
sequencing of data, provide a digital display to track 
successful counts of the Doppler frequency and, generally, 
control the flow of data between the signal processor and 
the Micro-VAX computer. Since a two component Laser Doppler 
Veloclmeter Is used In this experiment, the multiplexer 
sends four channels of data for each valid Doppler burst. 
The first data word contains the time Information, the 
second contains any status Information from the signal and 
the last two words contain the frequency data from each 
component. Due to the limited nature of this experiment, 
the first two words of data are not useful In further 
velocity calculations. Therefore, Just the frequencies from 
each particle caused Doppler burst are used to calculate the 
respective particle velocities. 

The Mlcro-VAX computer stores all acquired data in a 
buffer set up by acquisition programs. The next step after 
data acquisition Is to calculate the velocities. These 
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computer programs reduce the massive amounts of raw data 
obtained into meaningful velocity data that can be output in 
tabular form and graphical form. 

The following is a simplified path followed by the 
numerous computer programs to reduce the data. At least 1000 
samples from each component are required to obtain a statis- 
tically steady sample frequency. This frequency data must 
then be converted to the true Doppler frequency by the fol- 
lowing equation: 

f d = f s ‘ f B + f M 

where f d is the true frequency resulting from the particle 
traversing the probe volume, f g is the signal frequency 
measured by the counter, f 3 is the Bragg shift frequency 
of 40 MHz and f^ is the frequency of the downmixer. The 
downmlxer circuit permits decreasing of the effective freq- 
uency shift to values within the 1 to 10 MHz range. The true 
frequency data is then averaged over the samples taken to 
obtain a mean value and a sample standard deviation. Using 
this mean value, for each channel, all samples that are 
within 2 o-, two standard deviations, of the mean are saved 
and the rest are discarded. This procedure accounts for a- 
bout 5% reduction of data and will insure samples that are 
too low or high will not bias the final results. A new mean 
and standard deviation are calculated for each channel. Next 
the probe volume dimensions and parameters are calculated 
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with primary emphasis on the fringe spacing equation: 

« = X / [ 2 # sin(©/2>3 

where * Is the fringe spacing In microns, X Is the beam 
wavelength In microns and © Is the angle of Intersection 
of the beams. Knowing the fringe spacing, the velocity can 
be determined using the relation: 

u = f^ # 6 

where u Is the velocity of the particle In meters per second 
and f^ Is frequency of the particle measured by the count- 
ers In MHz. 

Once the velocity components of the particle are 
determined streamwlse and In the transverse directions, the 
computer stores the results In a data table for further 
analysis. The probe position Is also saved In this data 
file. The entire process is now repeated at different 
locations of Interest In the flow. 

C. LABORATORY APPARATUS 

Besides the Laser Doppler veloclmeter and schl leren 
systems described earlier, additional Instrumentation was 
used to obtain the wind tunnel pressure and temperature 
data. A 24 port scan 1 valve pressure transducer system was 
set up to measure the total and static pressures at various 
points Inside the operating tunnel. This scanlvalve system 
was driven by a digital to analog interface, built on the 
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premises, which allowed computer control of the "stepping" 
and "homing" of the scanlvalve. 

Using the Digital to Analog Interface on the MlcroVAX 
computer, these signals to step the scanlvalve system were 
transferred from a computer program to the physical unit for 
action. Once the stepper motor was sent the appropriate sig- 
nal to move, the Analog to Digital Interface board on the 
computer sampled the pressure transducer signals and loaded 
the data Into the computer for conversion Into pressure data 
and later analysis. 

The atmospheric pressure was recorded from a mercury 
barometer with micrometer scale mounted on the laboratory 
wall near the wind tunnel. The total temperature was record- 
ed on a mercury thermometer mounted on the barometer. Both 
of these measurements will soon be replaced with Paro- 
sclentlflc pressure sensors and thermocouples that will 
have a direct link with the computer. 

A TSI Six-Jet Atomizer, Model 9306, was used to provide 
particle seeding for the LDV system. 

Figure 2.9 shows the general laboratory setup and 
location of the data processing and ancillary equipment. 

D. DATA PROCESSING 

1 . Introduction 

The Laser Doppler veloclmeter yields a large amount 
of raw data during Its operation. This data must be stored 
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Figure 2.9 Laboratory Setup 



In the computer In a logical sequence In order to properly 
calculate the velocities and all other associated data. 

2. Computer Programs 



The foil owing 
tloning of each of the 
accomplish the task of 
the most efficient and 



discussion pertaining to the func- 
computer programs that was written to 
managing the vast amount of data In 
timely manner. 
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a. Master Program 



MASTER controls the overall flow of data during 
each of the subroutine calls, acquires all the data and 
reduces it to a usable form. 

First, the wind tunnel conditions are sampled by 
the use of various subroutines: DA11J1LAB, SETUP, GETDATA, 
DETACH and OUTPUTD. The traversing mechanism Is then posi- 
tioned with subroutine HWTRVT to move the probe volume of 
the LDV to the proper location in the flow. Next, the LDV 
data Is acquired with subroutine HWWA. The samples are 
taken, velocities calculated using subroutine VELOCITY and a 
statistical mean and standard deviation are calculated. The 
data Is then reduced with subroutine REDUCE_DATA to preclude 
the use of any velocities greater or less than two standard 
deviations away from the mean. Each of the two velocity com- 
ponents are then recorded along with the position in the 
flow. The velocity profiles are obtained from the data using 
post-acquisition processing programs. An option is given 
during the acquisition process to view a histogram of the 
resulting velocities to check on the distribution of data 
samples with subroutine HISTOGRAM. A sample histogram from a 
test run is shown in Figure 2.10. 
b. Scan 1 valve Programs 

First a signal is sent to the scanivalve to 
"home" to port 0 for the first pressure reading. This Is 
accomplished by using the subroutine DA11J1LAB. Next, a call 
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VERTICAL VELOCITY HISTOGRAM 




VELOCITIES (mps) 

Figure 2.10 Sample Histogram 



to subroutine SETUP is made to connect the computer's in- 
ternal timing clock and Analog to Digital interface board to 
program control. Also, the sample acquisition rate and num- 
ber of samples are interactively entered. 

Now, the computer is ready to acquire data from 
the scanivalve system. Since there are 24 ports to sample on 
the scanivalve, a loop is entered in which samples are 
taken, scanivalve stepped and data reduced during each loop. 
Subroutines GETDATA and DA11J1LAB acquire the data and step 
the scanivalve respectively. DA11J1LAB has an adjustable 
time delay to give time for the pressure in each port to 
equalize after scanivalve stepping. 
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After all 24 ports are sampled, the mean and 
standard deviation voltages are calculated for each port. 
The calibration pressure attached to port 0 and the atmos- 
pheric pressure at port 1 are now entered interactively in 
consistent units; inches of mercury in this case. With these 
two pressures entered, a calibration line is calculated and 
the remaining voltages from the scani valve system are con- 
verted to pressures. Using port 2, settling chamber stagna- 
tion pressure, and port 7, test section static pressure, a 
tunnel Mach number is calculated. See section III A for the 
equations used in these calculations. Next, the total 
temperature is entered and the free stream velocity is 
ca 1 cu 1 ated. 

For use in labelling a program generated plot- 
ting file for post-acquisition analysis, the angle of attack 
of the test section airfoil and current time are entered. 

The last part of the scanlvalve data acquisition 
programs close the loop on the analog to digital Interface 
by detaching the clock and driver. A diagnostic output file 
is also generated for post-acquisition analysis, 
c. LDV Data Acquisition Program 

Both the transmitting and receiving optics of 
the Laser Doppler Velocimeter are mounted on three axis 
motorized traverses that are slaved together to move the 
probe volume to any point of interest in the flow. Using a 
pre-existing software driver for the Taskmaster traverse 
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control unit, the optics could be moved at will from the 
computer workstation. 

Although a pre-existing software driver was 
available for the NASA-LDV Computer Interface, significant 
changes had to be made to the code for it to work on the 
present system. The main Laser Doppler Velocimetry 
acquisition subroutine "DAMUX" was responsible for loading 
the incoming data from the multiplexer into data arrays. 
Once the raw frequency data is in arrays, it has to be 
decoded to correspond to usable numbers. 

The time word output, as 14 bits, by the counter 
has a mantissa and exponent format: with a 10 bit mantissa 
and a 4 bit exponent. The computer reads it as a 16 bit word 
and adds two bits that are not used for data handling. These 
two bits are only used during direct computer links. This 16 
bit word is output as a complementary word that is reversed; 
that is, the zeros are output as ones and the ones as zeros 
then the entire word is reversed (i.e. the least significant 
bit which is on the right most side bit is moved to the left 
most bit in the word). First, the word is made into its com- 
plement then both the mantissa and the exponent are reversed 
prior to decoding. The exponent is preset from the counter 
front panel switch. The following is the standard 16 bit 
word composition and what actually results from the output: 
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Correct output . 



0111010001 0111 01 = 465x27 = 59,520 

Mantissa Exp NU 

Actual output 
01 0001 0111010001 
NU Exp Mantissa 

Translated by the computer 
0100010111 0100 01 = 279x24 = 4,464 

First step: Take Its complement 
10 1110 1000101110 
NU Exp Mantissa 



Next: reverse the order 
0111010001 0111 01 = 59,520 

Mantissa Exp NU 

where NU stands for Not Used. In order to fill the sixteen 
bit word, the 14th and 15th bit are added to control data 
flow but have no significance to the value of the data. From 
this last 16 bit word Is calculated the Doppler shift freq- 
uency. The output from the NASA-LDV Computer Interface Is 
actually representative of the time elapsed for the particle 
to traverse 8 fringes of the probe volume. The Doppler 
period Is expressed by the following equation: 
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t d = 1/32 CD 9 DQ...DQ) x 2^* ns 

where D 9 ..DQ represents the binary mantissa and EXP is 
the preselected exponent set on the counter. From the 
period, the Doppler frequency is calculated using the 
following equation: 

f d = 1000/t d MHz 



This Doppler frequency is then be used to calculate the 
velocity of the traversing particle. Since frequency 
shifting is used to discriminate the actual direction of the 
particle through the probe volume, the Doppler frequency 
must be adjusted to account for the Bragg and mixing 
frequencies. See equation on page 34 for this relationship. 
Once the actual frequency of the particle is determined, the 
fringe spacing is calculated by the first equation on page 
35 and the velocity is found using the second equation on 
the same page. 
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d. LDV Support Programs. 



Three subroutines are used to provide the 
supporting functions to the LDV data acquisition. These are 
VELOCITY, HISTOGRAM and REDUCE_DATA . 

VELOCITY is responsible for calculating the 
fringe spacing, extracting the particle frequency from the 
signal and then calculating the two velocities. 

HISTOGRAM takes the velocity data and plots 
a histogram of the percentage of velocities that fall within 
a given range. The x axis of the histogram is divided 
equally into 8 segments from the minimum to maximum velocity 
values. The option to plot a velocity histogram is given 
after each LDV acquisition to spot check the distribution of 
the data. A nearly Gaussian or normal distribution is 
expected if the system is operating properly. This will also 
show the effects of noise as evidenced by the broad band of 
the histogram values. 

REDUCE_DATA reduces the velocity data to 
only those values falling within two standard deviations of 
the mean. This results in a 4.5% rejection of data which 
will reduce the effect of high frequency noise on the 
counters and tend to reduce the effect of very large 
particles moving slowly through the probe volume. 
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III. DATA ACQUISITION 



A. TUNNEL CONDITIONS 

In order to obtain good correlation between the schlie- 
ren photographs and the Laser Doppler Veloclmeter obtained 
data, the operating tunnel conditions would have to be 
nearly Identical. Two sets of conditions were chosen: one 
corresponding to an approximate Mach number of 0.3 and the 
other 0.5. Once the tunnel was running at operating 
conditions, the computer program MASTER calculates the test 
section Mach number and free stream velocity according to 
the following relationships. Assuming isentroplc flow in the 
test section, the section Mach number, M, is calculated from 
the stagnation pressure, P Q , and the static pressure, P, 
using the following relationship: 

M * (([P 0 /p] ( y” 1/y> -l)x<2/[y-l])) 1/2 

where y = 1.4. The resulting calculated Mach numbers for 

the two conditions are as follows: 0.295 and 0.450. From 

these Mach numbers, the corresponding free stream velocities 
are calculated. 

V = M # a 

where a is the sonic velocity. 

a - <yRT> 1/2 

where R is the universal gas constant and T is the static 
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temperature In the test section. The static temperature can 



be calculated from the total temperature, T Q , measured on 
the wall near the tunnel by the next relationship: 

T = Tq/U + [ (y-l >/23#M 2 ) 

B. SCHLIEREN PHOTOGRAPHS 

1 . Introduct 1 on 

Once the wind tunnel was running at the first of 
two operating conditions, manometer board readings were 
recorded and Mach number calculated using the first 
equation on the prelous page. The schlleren photographs were 
then taken at four different angles of attack and three 
different strobe Intensities to obtain the optimal 
photographs . 

2. Test Conditions and Procedure 

Using the procedures detailed In section II A, the 
schlleren photographs were taken at the settings and 
conditions listed In Table I, below. 

Upon completion of the wind tunnel runs, the 35mm 
film was developed and printed. From these schlleren 
photographs, the conditions resulting In the best depiction 
of the wake profiles were chosen to represent each Mach 
number and angle of attack. These are annotated with an 
asterisk, #, In Table I, below. 

Using these photographs, a short HP-41 CX calculator 
program was written to relate photographic dimensions to 
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TABLE I. 


SCHLIEREN PHOTOGRAPH PARAMETERS 




n 


a (deg) 


Strobe 

Intensl tv (Hz 


# 


0.30 


0.0 


20 




0.30 


0.0 


70 




0.30 


0.0 


150 




0.30 


7.7 


20 




0.30 


7.7 


70 


# 


0.30 


7.7 


150 


# 


0.30 


15.5 


20 




0.30 


15.5 


70 




0.30 


15.5 


150 




0.30 


23.6 


20 




0.30 


23.6 


70 




0.30 


23.6 


150 


# 


0.45 


0.0 


20 




0.45 


0.0 


70 




0.45 


0.0 


150 


# 


0.45 


7.7 


20 




0.45 


7.7 


70 




0.45 


7.7 


150 




0.45 


15.5 


20 


# 


0.45 


15.5 


70 




0.45 


15.5 


150 



real test section dimensions. Once the real dimensions of 
the wake were known, the Taskmaster traverse mechanism was 
calibrated. This yielded a 1:1 relationship in the chordwise 
direction and a 1:4 relationship in the vertical direction. 
This means that the Taskmaster digital control panel reads 
true values for the X direction, chordwise, and four times 
the real value in the Z direction, vertical. Arbitrarily 
choosing the trailing edge of the airfoil at 0° a as an 
initialization point for the Taskmaster, a table of LDV 
sampling points was made for the first three angles of 
attack. See Table II, below, for a sample of the Taskmaster 
table at a = 7.7° 
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Table II. TASKMASTER LDV SAMPLING LOCATIONS 
Points Alpha = 7.7 deg 



Sampled £ 


l 


~> 
































20 2.00 


1 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 


2 


1.5 


1 


.5 


0 


-.5 -1 


-1.5 -2 -2.5 


20 1.50 


7 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 


2 


1.5 


1 


.5 


0 


-.5 -1 


-1.5 -2 -2.5 


20 1.00 


7 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 


2 


1.5 


1 


.5 


0 


-.5 -1 


-1.5 -2 -2.5 


20 .50 


7 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 


2 


1.5 


1 


.5 


0 


-.5 -1 


-1.5 -2 -2.5 


20 .00 


7 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 


2 


1.5 


1 


.5 


0 


-.5 -1 


-1.5 -2 -2.5 


19 -.50 


7 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 


2 


1.5 


1 


.5 


0 


-.5 


-1.5 -2 -2.5 


17 -1.00 


7 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 


2 


1.5 


1 


.5 






-1.5 -2 -2.5 


17 -1.50 


7 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 


2 


1.5 


1 


.5 






-1.5 -2 -2.5 


11 -2.25 


7 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 














-2.5 


12 -3.00 


7 


6.5 


6 


5.5 


5 


4.5 


4 


3.5 


3 


2.5 


2 












-2.5 



176 (Total locations sampled) 



C. LASER DOPPLER VELOCIMETRY DATA ACQUISITION 
1 . Calibration of the LDV System 

The relationship between the Doppler frequency and 
the velocity of the'partlcle traversing the probe volume Is 
given by the second equation on page 35. The particle 
velocity is directly proportional to the Doppler frequency 
by a factor of the fringe spacing. By the first equation on 
page 35, fringe spacing is only a function of the 
transmitting lens angle and wavelength of the light. Since 
the wavelength of the light is accurately known to within 
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0.05 microns, or 0.01%, an accurate lens angle Is all that 
is needed for calibration [Ref. 4], 



The following relationship was used to calculate 
the lens angle: 

e = 2 # tan -1 <E # d / ( 2 # f)> 

where E is the beam expansion ratio, 1.0 in this case: d is 

the incident beam spacing, preset at 13 mm with TSI Model 

9113-13; and f is the wavelength of the light. 

The lens angle in this case was not accurate 

because the focal length of the focusing lens was not 

known. Therefore, calibration was necessary to Insure 

accurate velocity measurements were obtained. 

The calibration method chosen was to operate the 
wind tunnel at known conditions, calculate the free stream 
velocity from the scan! valve pressure data and then move 
the probe volume via the Taskmaster to coincide with a 

position in which the free stream velocity could be 

measured. From this procedure, a constant could be applied 
to the velocity equation as follows: 

u = f^ # 6 # (Calibration Constant) 

This calibration constant will compensate for the estimated 
focusing lens focal length. 

2. Test. Conditions and Procedure 

Once calibration is completed, the tunnel is 

brought to the test conditions, the Taskmaster is moved to 
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the first sample point and the MASTER computer program is 
run with the adjusted fringe spacing to reflect the 
cal ibrat i on . 

At each sampling point, 1000 samples will be taken 
to Insure statistical accuracy of the mean flow calcula- 
tions. The scheduled sampling points can be seen in Table 
II, page 48. After each Taskmaster traverse, the oscillo- 
scope is checked to Insure no drastic loss of the signal 
occurred during the movement. If the Doppler signal de- 
creases, slight adjustments to the receiving optics traverse 
may be needed to regain the signal strength. To optimize the 
data received, slight adjustments may have to be made, also, 
to the selected mixing frequency, mixing frequency ampli- 
tude, counter threshold setting or high and low pass filter 
sett i ngs. 

After Table II data points are sampled, the airfoil 
angle of attack or test Mach number is changed and the above 
procedure is repeated. 
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IV. RESULTS AND DISCUSSION 



A. SCHLIEREN 

The first thesis goal was attained early enough to 
provide qualitative Information concerning the airflow 
around and In the wake of the airfoil. Schlleren 35mm 
photographs were produced for both selected mach numbers, 
0.3 and 0.45, and for four different angles of attack; 0°, 
7.7°, 15.5° and 23.6°. The best photographs for each 
Mach number and the first three angles of attack are dis- 
played In Appendix A, pages 55-57. From these photographs, 
LDV sampling points were chosen as seen from Table II on 
page 48. 

Three turbulent areas are evident on all the schlleren 
photographs. As seen In Figure A. 4 on page 56, the airflow 
Is slightly disturbed upstream of the airfoil due to the 
presence of a pitot-static tube extending Into the airflow 
Just forward of the airfoil. This tube has been subsequently 
removed to clean up the flow. At the top of the same photo- 
graph, a wedge shaped disturbance can be seen. This and the 
lower floor turbulent layer are results of the transonic 
wall slottings found on both the roof and the floor of the 
test section. Although a concern for later unsteady flow 
studies, these two areas do not Interfere with the wake 
profiles resulting from a steady flow analysis. 
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B. LASER DOPPLER VELOCIMETER 



Due to equipment problems and time constraints, the LDV 
data acquisition part of this thesis was not able to be 
accomplished. All systems are now operational, to include 
computer programs, however, there is Inadequate time to 
acquire data and analyze the velocity profiles obtained. 

Calibration runs were made with a limited degree of 
success when the wind tunnel was run under Mach 0.3 and the 
Laser was not powered over 2 watts. The problem with the 
Laser is easily repaired, however, time consuming. The 
upper limit on the wind tunnel Mach number is caused by the 
evaporation of the seeding particles in high speed flow. 
This, too, can be easily overcome by the proper selection 
of seeding particles and related atomizer. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The thesis provided invaluable assistance to the dynamic 
stall project at the Fluid Mechanics Laboratory. The schlie- 
ren photography was able to provide the qualitative data on 
the steady state flow field. With this accomplished, the 
schlieren system in this laboratory was able to be easily 
tailored to the next step: unsteady flow analysis. 

The LDV work provided the laboratory with the essential 
computer programs, hardware interfaces and drivers, optics 
systems and auxiliary required equipment to enable the pro- 
ject to proceed with very little additional work to achieve 
the steady state velocity profiles. The ground work has been 
accomplished to enable the project to continue forward with 
the systems on hand to study the effects of dynamic stall in 
both the steady and unsteady environments. 

B. RECOMMENDATIONS 

The schlieren system is fully functional, however, the 
Laser Doppler Velocimeter system still requires a little 
additional effort to make it productive. The Laser doughnut 
condition needs to be addressed to enable higher powers to 
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be utilized. The types of seeding ayailable to the system 
needs to be addressed to insure operability of the tunnel 
and data acquisition systems above Mach 0.3. With these 
minor fixes and a little fine adjustments of the signal 
processing equipment, the LDV is ready for action. 
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APPENDIX A 



SCHLIEREN PHOTOGRAPHS 




Figure A.l Mach Number = 0.3, a = 0° 




Figure A. 2 Mach Number = 0.45, a = 0° 
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Figure A. 6 Mach Number = 0.45, a 



15.5° 
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APPENDIX B: COMPUTER PROGRAMS 



PROGRAM MASTER 



Master Progrom written by Maj Richard R Ryles, NPS Student 28 Apr 87 

Last Update 6 Sep 87 

This program will gather data from the scanivalve system, move the 
Taskmaster and acquire LDV velocity data. 

This progrom should be used with subroutine SETUP. GETDATA. DETACH. OUTPUTD 
HWTRVT. HWWA. HISTOGRAM ond VELOCITY. 



INCLUDE ’ SYS$LIBRARY : LIOSET . FOR ’ 



Declare variables 



INTEGER DEV_ID.CLX_ID 

INTEGER NCHANL.HISTO.BLUE_GREEN 

INTEGER NPT, NPTS.NSAiJ»LES. COUNT 

INTEGER I .J.NSTEPS. IER. ANSWER. END. TALLY. TIME 

REAL VOLTAGES(5000).U_U(200).V_V(20O) , U_V(200) .DELX.DELZ 

REAL'S SUM . AMEAN . X . XX . Y . YY . STD.DEV . STD_DEV_V . STD_DEV_U . SUM_V . SUM.U 

REAL'S SUM_VERT,SUM_HORIZ 

REAL RAT E . SCAN I _VOLT (0:23.2) . SCANI .PRESS (0: 23) . VERT_VEL(8200) 
REAL CALIB.PRESS. SLOPE. INTERCEPT. HORIZ_VEL(8200).FINAL_V.FINAL_U 
REAL X.POSITION, Z_POSITION.MEAN(2) ,M. FREE_STREAM, TO, ATMOS_PRESS 
REAL X_LOC(200) ,Z_LOC(200) ,U(200) ,V(200) .ALPHA. ZERO. CORR.FACTOR 
REAL'S UU.W.UV 
CHARACTER' 1 ANS.VOLTJM6 
LOGICAL HOME 



Scanivalve dato acquisition 



I 



OP EN(UNIT-99. FILE-’ RESULTS. DAT ’.STATUS- ’NEW’) 

TALLY - 0 
HOME - .TRUE. 

NSTEPS - 0 

CALL DA1 1 J 1 LAB (HOME. MST EPS. IER) 

IF (IER.NE. 1) GOTO 17* 

ANSWER - 0 

WRITE ('.') ’Skip to Taskmaster? (y/n) [DEF - No]’ 

READ (*.1) ANS 

IF (ANS.EQ. ’Y’ .OR.ANS.EQ. *y’) GOTO 90 

WRITE ('.') ’Print th# scanivalve diagnostic data to a file? (Y/N)’, 
k ’ (DEF - Mo]’ 

READ ('.1) VOLT.ANS 

CALL SETUP (CLK_ID.DEV_ID.RATE.NCHAML.NPT.NPTS) 

0 CONTINUE 

HOME - .FALSE. 

NSTEPS - 1 
DO 50 J - 0. 23 
WRITE ('.2) J 

IF (J.NE.0) CALL DA1 1 J 1 LAB (HOME. MST EPS. IER) 

IF ( IER.NE. 1 ) GOTO 150 
CONTINUE 

CALL GETDATA (DEV_ID .NPT . NPTS .VOLTAGES) 



Calculate neon and stondord deviation 



SUM - 0.0 

DO 30 1-2. NPT 

SUM - SUM + VOLTAGES(I) 
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30 CONTINUE 

AMEAN - SUM / FLOAT (NPT-1 ) 

SUM - 6.M 
DO 46 I -2, NPT 

X - VOLTAGES! I) - AMEAN 
Y - X • X • (1 .e/FLOAT (NPT-1)) 

46 SUM - Y + SUM 

STD.DEV - SORT (SUM) 

IF (VOLT.ANS.EQ. ’y’ .OR.VOLT.ANS. EQ. ’Y’) CALL OUTPUTO (RATE. NPT. 
k NCHANL. AMEAN. STD.DEV, J. VOLTAGES) 

SCANI_VOLT(J.1) - AMEAN 
56 SC AN I _ VOLT (J. 2) - STD.DEV 

CALL DETACH (CLK_ID.DEV.ID) 

IF(VOLT_ANS.EQ. ’Y’ .OR . VOLT.ANS. EQ. ’y’ ) THEN 

OPEN ( UN I T- 1 5 . F I LE- * SCAN I.VOLT . DAT ’ . STATUS- * NEW * ) 

DO 60 I - 0. 23 

60 WRITE (15.3) I. SCAN I.VOLT ( I . 1 ) . SCAN I.VOLT (1.2) 

END IF 



C Calculate scanivalve calibration line 
C and convert all voltages to pressure 



WRITE (*.*) ’Enter Scanivolve Calibration Pressure:' 

READ (*,*) CALIB.PRESS 

WRITE( *,*)’ Enter Calibration Pressure Correction Foctor:’ 

READ(* , • )CORR_FACTOR 

CALIB.PRESS - CALIB.PRESS + CORR_ FACTOR 
WRITE(* ,*)’ Enter atmospheric pressure:’ 

READ(* , • ) ATUOS.PRESS 

SLOPE - ( ATMOS.PRESS— CALI B.PRESS)/ (SCAN I.VOLT ( 1 , 1 )-SCANI_VOLT (0 . 1 ) ) 

INTERCEPT - CALIB.PRESS - SCANI_VOLT(6. 1 ). SLOPE 

WRITE(99, •) ’CALIB.PRESS ’ .CALIB.PRESS 

WRITE(99, •) ’CORR.F ACTOR ’ , CORR.FACTOR 

WRITE(99, • ) ’ ATMOS.PRESS ’ , ATMDS.PRESS 

OPEN ( UN I T- 1 6 . F I LE- * SCAN I PORT . DAT ’ . STATUS- ’NEW’) 

DO 76 I - 0. 23 

SCAN I .PRESS (I) - SCAN I .VOLT ( I , 1 ) • SLOPE + INTERCEPT 
70 WRITE(16, »)I , SCANI.PRESS(I) 

CLOSE( 16) 

IF (VOLT.ANS.EQ. *Y’ .OR.VOLTJtNS.EQ. ’/’) THEN 
DO 86 I - 0. 23 

80 WRITE (15.4) I .SCAN I .PRESS ( I ) .SCAN I .VOLT ( 1 . 2) *SLOPE+INTERCEPT 

END IF 
CL0SEO5) 

WRITE!*.*)’ ’ 

WRITE!* . • ) 'Equat ion of the colibrotion line -‘ 

WRITE(*.*)‘Y -’.slope.’ X +'. intercept 
WRITE(99.*)’Y -’.slope.’ X ♦’ . i ntercept 
WRITE(* , •) ’ ’ 



C Calculate tunnel Wloch and velocity 



I F(SCANI_PRESS(7) . EQ. 0 . 0) SCANI_PRESS(7) - 1.0 
IF(SCANI_PRESS(2) . LT . SCANI_PRESS(7) ) GOTO 85 
M - SQRT( ( (SCANI_PRESS(2)/SCANI_PRESS(7) ) •• (2 ./7. )— 1 . )*5. ) 

85 WRITE(* , • ) ’Tunne I Mach Number -’ ,M 
WRITE(99, •) ’Tunnel Moch Number -’ ,M 
WRITE( * . • ) ’ * 

WRITE( *,*)' Enter Total Temperoture (C)’ 

READ( • , • ) TO 

WRITE(99.*) 'TOTAL TEMP’, TO 
WRITE( * . • ) ’ ’ 

TO - TO + 273.15 

FREE.STREAM - M • SQRT(401 ,52*TO/(1 .+(M**2/5. ))) 

WRITE!* ,*) ’Free Stream Velocity (mps) -’.FREE.STREAM 
WRITE(99.*)’Free Stream Velocity (mps) -’.FREE.STREAM 
IF(FREE.STREAM.EQ.0.0) FREE.STREAM - 1.E-10 
WRI TE( * . • ) ’ ' 
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00000 — 000 OOO — <0000 



WRITE( •,•)’ Enter angle of attack (deg):' 

R£AO( • . • ) ALPHA 
WRITE(99 , •) ’ALPHA’ . ALPHA 

WRITE(» ,•) ’Enter Current Time: (4 digit integer)’ 
Reod(*,*) TIME 
WRITE(99. •) ’TIME’ .TIME 

WRITE (•.•) ’Skip to end? (y/n) [DEF - Mo)’ 

READ (*,1) ANS 

IF (ANS. EG. ’y* .OR. ANS. EG. ’Y’) GOTO 100 



Taskmaster mover and LDV data acquisition 



I 



WRITE (*.*) ’Which color is connected to the left (#3) channel?’, 
t ’ (B I ue-1 .Green-©) ’ 

READ (*.9) BLUE_GREEN 
WRITE(99 , • ) ’BLUE_GREEN’ ,BLUE_GREEN 
IF(FREE_STREAM. EO. 0. 0) FREE.STREAM - 1 . E-0 
OPEN (UNIT-25. FILE- ’FINAL. DAT’ . STATUS- ’ H EW ’ ) 

WRITE (25, •) ’ Taskmaster (Final. Dot)’ 

WRITE (25, •) ’ pos (in) velocity/std dev velocity/std dev’ 

t«ITE (25.16) 

WRITE (•.•) ’Skip to LDV? (y/n) [Def - NO]’ 

READ ( • . 1 ) ANS 

IF ( ANS . EO . ’ Y ’ . OR . ANS . EQ . ’ y ’ ) GOTO 110 
IF(TALLY.EQ.e) THEN 

WRITE( •,•)’ Enter X-Position for calibration’ 

READ( • . • ) X_POS I T I ON 

WRITE(* ,•)’ Enter Z-Position for calibration’ 

READ( • , • ) Z.POSITION 
ENO I F 

CALL HWTRVT (DELX.DELZ) 

WRITE (•,•) ’Move the Taskmaster again prior to LDV Acquisition? (Y/N)’ 
READ (*.1) ANS 



Record position of LDV at all times 



I F ( ANS . EO . ’ Y * . OR . ANS . EQ . ’ y ’ ) THEN 
X POSITION - X.POSITION + DELX 

Z.POSITION - Z.POSITION + DELZ 
GOTO 100 

EX) IF 



Acqu i re LDV data 

10 CALL HWWA (NSAAPLES. ANSWER. BLUE.GREEN) 



Increment position and record it 
Read velocity data, sum it and find 
standard deviation 



X POSITION - X_POS I T I ON + DELX 
Z POSITION - Z.POSITION + DELZ 
SUM. VERT - 0.D0 
SUM HORIZ - O.D0 

OPEN (UNIT-35. FILE-’VELOCITY. DAT’ . STATUS- ’ OLD ’ ) 

DO 120 I - 1. NSAAPLES 

IF (BLUE GREEN. EO.0) READ(35.7) VERT_VEL( I ) ,HORIZ_VEL( I ) 
IF (BLUE.GREEN . EQ. 1 ) READ(35,7) HORIZ_VEL( I ) . VERT_VEL( I ) 
SUM VERT - SUM.VERT + VERT.VEL(I) 

120 SUM.HORIZ - SUM. HORIZ + HORIZ.VEL(I) 

CLOSE(35) 

FINAL V - SUM.VERT / FLOAT (NSAAFLES) 

FINAL.U - SUM.HORIZ / FLOAT (NSAAFLES) 

SUM.V - 0.D0 
SUM.U - 0.D0 
UU - 0 D0 
W - 0 00 
UV - 0 D0 



60 



DO 13# I • 1 , NSAO>LES 

X - VERT VEL(I) - FINAL.V 

XX - HORI Z.VEL ( I ) - FINAI U 

Y - X • X • ( 1 . 6/FLOAT (MSAJkFLES-1 ) ) 

YY - XX • XX • (1 .0/FLOAT(NSAJ*>LES-1)) 

UU - UU + ((HORIZ_VEL( I )**2)-(FINAL_U»»2)) 

W - W + ( (VERT.VEL ( I )• •2)-(FINAL_V»»2) ) 

UV - UV ♦ ((HORIZ_VEL(I)»VERT_VEL(I))-(FINAL_U»FINAL_V)) 
SUM_V - Y + SUM.V 
130 SUM.U - YY + SUM.U 

UU - UU / FLOAT (NSAfcPLES) 

W - W / FLOAT (NSAMPLES) 

UV - UV / FLOAT (NSAfcPLES) 

STD_DEV_V - SORT (SUM.V) 

STD_DEV_U - SORT (SUM.U) 

WRITE(* , • ) * ’ 

WRITE (• , • VVeloci t ies ond std devs -’ 

WRITE (• , •) 'U dir’ , FINAL.U , STD.DEV.U 
WRITE l • , •) ’V d i r ’ , FINAL.V, STD.DEV.V 
WRITE (99. •) 'Veloei t ies ond std devs -’ 

WRITE (99, •) ’U di r’ , FINAL.U, STD.DEV.U 
WRITE i 99, •) ’V dir’, FI NAL_V , STD_DEV_V 
WRITE (•.•)’ ’ 

WRITE(* , • ) ’ Fi rst cut RMS volues’ 

WRITEi* , •) ’ U-RMS , V-RMS ’ ,UU, W 

WRITEi*, •) ’UV-RMS’ ,UV 

WRITEi 99, •) 'Fi rst cut RMS volues’ 

WRITE(99, •)’ U-RMS, V-RMS ’ ,UU,W 
WRITE(99,*)’ UV-RMS’ ,UV 



C Eliminote doto outside two sigmo from meon 



IF (BLUE.GREEN.EQ.1) THEN 
MEAN(1) - FINAL_U 
MEAN(2) - FINAI V 

CALL REDUCE.DATA ( NSAJyPLES , COUNT , END , MEAN , STD.DEV.U , STD.DEV.V, 
t BLUE.GREEN) 

FINAI U - MEAN( 1 ) 

FINAL.V - MEAN(2) 

WRITE (99, •) ’Reduced velocities ond std devs «’ 

WRITE (99 , • l ’U dir’ ,FINAL_U.STD_DEV_U 
WRITE (99 , •) ’ V di r’ ,FINAL_V.STD_DEV_V 

ELSE 

MEAN(2) - FINAL.U 
MEAN( 1 ) - FINAL.V 

CALL REDUCE.OATA (NSAMPLES, COUNT. END. MEAN. STD_DEV_V,STD_DEV_U, 
It BLUE.GREEN) 

FINAL.U - MEAN ( 2 1 
FINAL.V - MEAN(1) 

WRITE (99,*) 'Reduced velocities ond std devs «’ 

WRITE (99, *1 ’U dir’, FINAL_U,STD_DEV_U 
WRITE (99,*)’V di r ’ . FINAL.V, STD.DEV.V 

ENOIF 



C Colculote RMS volues 



OPEN(UNI T-43 , FI LE-’VELOCI TY . DAT ’ . STATUS- ’OLD ’ ) 

UU - 0. DO 
W - 0.DO 
UV - 0.D0 
00 210 I - 1. END 

IF (BLUE.GREEN.EO.0) READ(4J,7) VERT.VEL ( I ) , HOR I Z.VEL ( I ) 

IF (BLUE.GREEN.EQ.1) READ(43,7) HOR I Z.VEL ( I ) , VERT.VEL ( I ) 

UU - UU -f ((HORIZ_VEL( I )**2)-( FINAL.U* *2)) 

W - W + ( ( VERT.VEL ( I ) • *2)-( FINAL_V**2) ) 

210 UV - UV -f ( (HORI Z.VEL ( I ) ‘VERT.VEL ( I ) )-( FINAL.U* FINAL.V) ) 

CL0SE(43) 

UU - UU / FLOAT ( EM) ) 
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W - W / FLOAT (DO) 

UV - UV / FLOAT (DO) 

WRITE(* , •) ' * 

WRITE! • , •) ‘Reduced RMS values' 

WRITE!*,*) ’ U-RMS , V— RMS ’ .UU.W 
WRITE!*,*) ‘UV-RMS* .UV 
WRITE! 99 . •) 1 U-RMS . V-RMS ' .UU.W 
WRITE(99.*)'UV-RMS‘ .UV 

WRITE (25.8) X.POS IT ION. Z.POSIT ION. FINALJJ.STD.OEV.U. FINALLY. STD DEV V. 
k UU.W.UV 

TALLY - TALLY + 1 



Options to vis* histogram ond toks doto ogoin 



WRITE(* . •) ’ * 

WRITE (*,*) ‘Do you wont o Histogram? (1/0)’ 

READ (•.•) HISTO 

IF (HISTO. EQ.1) CALL HISTOGRAM (COUNT .END. BLUE GREEN) 
WRITE (*.6) 

READ (*.1) ANS 

IF (ANS.EQ. ‘n'.OR.ANS.EQ. ‘N’) GOTO 149 
GOTO 100 
40 CLOSE (25) 

CL0SE(99) 



Convsrt output dato into propsr formot for 
"PLOTCV" Vector Plot of velocity profiles 



OP EN(UNIT-45. FILE- ‘FINAL. DAT* .STATUS- ‘OLD’) 

READ! 45,14) 

OPEN ( UN I T-46 . F I LE- ’ GRA . DAT ‘ , STATUS- ’ NEW ‘ ) 

WRITE(46. 13) U,FREE_STREAM. ALPHA, TIME, 11 .16.10 
ZERO - 0.0 
DO 150 1-1. TALLY 

READ (45, 1 1 )X_LOC( I ) , Z_LOC( I ),U(I),V(I) ,U_U( I ) , V_V( I ) ,U_V( I ) 

X_LOC( I )-X_LOC( I )+3 . 0 
Z_LOC( I )-(Z_LOC( I )+2 . 5 ) * . 24 
U(I) - U(I) / FREE.STREAM 
V(I) - V(I) / FREE.STREAM 
U_U( I ) - U_U(I) / FREE_STREAM 
V_V(I) - V_V(I) / FREE_STREAM 
U_V(I) - U_V(I) / FREE_STREAM 

WRITE(46 . 1 2 ) I . X_LOC( I ) . Z_LOC( I ) . U( I ) . V( I ) . ZERO.U_U( I ) . V_V( I ) . ZERO. 
k U_V( I ) . ZERO 

CLOSE(45) 

CLOSE(46) 

WRITE ( • . 5) 



Formot statements 



I 



FORMAT (A1 ) 

FORMAT (IX. ‘Sconivol vs f M2) 

FORMAT (IX. 'Port 12 . 3X , ‘Average Voltage -'.F12.4.3X, 
k ‘Standord Deviation -’.F12.4) 

FORMAT (IX. ‘Port f ‘ . 12 . 3X. ‘Pressure -\F12.4,* (psi)'.3X, 
k ‘Standard Deviation -'.F12.4,* (psi)’) 

FORMAT (//IX. ‘This is the ENOI ’//) 

FORMAT (/IX, ‘Move taskmaster ogoin and toke dato? (y/n) [DEF - Y]’/) 
FORMAT (1X.2F11.4) 

FORMAT (1X,2F7.2,F10.2,F8.2,F10.2,F8.2,1 P3E12 .4) 

FORMAT (II) 

11 FORMAT (IX, 2F7 .2,F10.2,8X,F10.2,8X, 1P3E12 . 4) 

12 FORMAT (1X.I3.1P10E12.4) 

13 FORMAT ( 1 X , ' GRA . DAT , MACH MO - *.F4.2.‘. FREE STREAM VELOCITY - ‘.F5.1. 

k *. ALPHA - ' ,F4. 1 , /IX. ‘September 9, 1987* . 4* . i 4 . * hrs. ‘ 

k //4X.’H>ROF NSAU>L NCOEF ‘/7X . 12 . 6X . 12 . 6X . I2/1X . * ID X-POS* . 

k ' Z-POS U V W Ut2 ' , 

k Vt2 Wt2 UV UW‘ ) 
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16 



F0RU4T //) 

FORMAT (0X, ,0X, * 2 * . 12X, *u* , 17X, *v* . 1 IX. *ut2* ,9X, * vt2’ . 10X, *uv* ) 
STOP 



C 

C Sconlvolve Error Statements | 

C 

170 IF(HCME) THEN 

•(RITE (»,•) ’Problom in honing Scanivalves' 

ELSE 

WRITE (• , • ) ’Probl on in stopping Scanivalves* 

ENOIF 

IF (IER.EQ.2) WRITE (•.•)’Dota tronsfor failod in J-8oord’ 

IF (IER.EQ.3) WRITE (•,•)* Fa i I lire in detaching J-Boord* 

IF (HOME) GOTO 10 

GOTO 20 

END 
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SUBROUTINE SETUP(CLK_ID.DEV_ID.RATE.MCHANL.hPT.NPTS) 



C Subroutine written by Maj Richard R Ryles, KPS Student, 13 May 87 
C Updoted 15 Jul 87 

C 

C This subroutine is intended to ottoch the internol clock ond A-0 boord 
C to the mo i n program ond allow them to be used for doto acquisition. 

C 

C This subroutine should be used with subroutines GETDATA, DETACH 

C 

INCLUOE •SYSSLIBRARY: LIOSET.FOR* 

C 

C Declore locol voriobles 



CHARACTER* 4 CLK_DEV 
CHARACTER* 1 ANS 






INTEGER SYS.STAT 


IStotus ret 


by LIO calls 


INTEGER DEV_ID,CLK_ID 


! LIO dev i ce 


ID 


INTEGER NCHAN 
INTEGER NPT, NPTS 
REAL RATE. TIME 


! Chonne 1 f 




Input desired clock ond 


sample parameters 


1 



WRITE( • . 1 ) 

READ( • . 2 )NCHANl 
10 WR I TE( • , 3) 

READ( • , • )RATE 
I F(RATE. GT . 600 . 00) THEN 
WRITE(*,4) 

GOTO 10 

ELSE 

IF(RATE.LT.0.01)WRITE(* .5) 

I F(RATE. LT . 0 . 01 )GOTO 10 

END IF 

20 WR I TE( • , 6) 

READ( • , 7) NPT 
I F(NPT . GT . 5000)THEN 
WR I TE( • ,8) 

GOTO 20 

ELSE 

IF(NPT.LE.0)WRITE(*,9) 

I F(NPT . LE . 0)GOTO 20 

END IF 

NPT - NPT + 1 

TIME - FLOAT ( NPT- 1)/RATE 

WR I TE( • . 1 1 )TIME. TIME/60. 

NPTS - 2 • NPT 
CLK DEV - ’ KZA0 ' 

WRITE(*. 12) 

READ( • . 1 3)ANS 

IF ( ANS . EO . ’ y ’ . OR . ANS . EO . ’ Y ' ) THEN 
WRITE(* . 14) 

READ( • . 1 5)CLK_D£V 

END IF 

c 

C Attach the KWV11— C Clock ond set up for mopped I/O 

C 

SYS.STAT - LIOSATTACH (CLK_ID, CLK_DEV. LIO$K_MAP) 

I F ( . NOT . SYS STAT) CALL LIB$SIGNAL(XVAL(SYS_STAT)) 
SYS.STAT - LIOSSET.R (CLK_ID. LIO$K_CLK_RATE. 1 .RATE) 

I F( . NOT . SYS.STAT) CALL LIB$SIGNAL(XVAL(SYS_STAT)) 
SYS.STAT - LIOSSET.I (CLK.ID. LlO$K_TRIG. 1. LI0$K_II*4EDIATE) 
I F ( . NOT . SYS.STAT) CALL L I B$S IGNAL(XVAL (SYS.STAT ) ) 

C 

C Attoch the ADV11-D ond set up for mopped I/O 

c 

SYS.STAT - LIOSATTACH (DEV_ID. ’ AZA0 * . LIO$K_UAP) 



64 



ooo oooooo 



r - • 



I F ( , NOT . SYS.STAT) CALL L I B$S I CNAL(%VAL( SYS.STAT) ) 



Set op the ADV for synchronous transfer | 

Sat chonnel (nchonl) 

Coin of 1 

Trigger on LIO$READ and fill buffer as fast as possible 



SYS.STAT- LIOSSET.I (DEV.ID. L I OSK.AD.CHAN . 1 , NCHANL, ) 

I F ( . NOT . SYS.STAT) CALL L I B$S I GNAL(XVAL( SYS.STAT ) ) 

SYS.STAT - LI OSSET. I (DEV.ID. LlO$K_AD_GAIN , 1. 1.) 

I F ( . NOT . SYS.STAT) CALL LIB$SIGNAL(XVAL(SYS_STAT) ) 
SYS.STAT - L I OSSET. 1 (DEV.ID, LIOSK.TR IG . 2 . LIOSK.CLK.POINT.CLK.ID) 
I F ( . NOT . SYS.STAT) CALL LIB$SIGNAL(XVAL(SYS_STAT) ) 



Format Statements 



1 FORMAT ( 1 x En t e r channel number to sample: (0, 1,8,9)’) 

2 FORMAT (12) 

3 F0RMAT(/1x, ’ Input desired Clock Rote (600 Hz - 0.01 Hz):’/) 

4 FORMAT( 1 x , ’ Fo I I ow i ng directions! Input too h i gh . ’//20x . ’ Try again'/) 

5 FORMAT( 1 x , ’ Too low! Lots of time on your hands? ’//20x ,' T ry ogoin’/) 

6 F0RMAT(/1 x . ’ No . of doto points: (any number less than 5,001)'/) 

7 F0RMAT(I4) 

8 FORMAT (/I x ,’ Hey ! Follow directions! Too high. Try ogoin!’/) 

9 F0RMAT(/1 x . ’ Hey ! Lost input is nonsense! Try ogoin!'/) 

11 FORMAT (/l x ,’ The opproximote run time for each port will be * . f 6 . 1 , 

k ’ seconds (or ’,f5.2.’ minutes)’/) 

12 FORMAT (/I x ,’ Do you wont to chonge the Clock Device Nome from [KZA0]?’, 

4 • (y/n) [DEF - no]’/) 

13 FORMAT ( A 1 ) 

14 FORMAT( lx, ’Enter Clock Device Nome: (KZA0 , KZB0) ’ ) 

15 F0RMAT(A4) 

C 

RETURN 

END 



1 
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ooo oooooo ooo ooo ooooooooooo 



SUBROUTINE G ETDATA(DEV_ ID, NPT ,NPTS, VOLTAGES) 



Subroutine written by Moj Richard R Ryles. NPS Student, 28 Apr 87 I 

Update 15 Jul 87 j 

I 

This subroutine it intended to ocquire the dato front the A-0 boord 
with the specified parameters from subroutine SETUP. It olso will | 

translate the A-0 data into voltages to be converted to physical j 

quantities in the main program. 

I 

This subroutine should be used with subroutines SETUP , DETACH . OUTPUTD j 



INCLUDE * SYS$ LIBRARY : L IOSET . FOR * 



Declare local variables 



i 



INTEGER SYS.STAT 
INTEGER DEV_ ID 
INTEGER DATA.LENGTH 
INTEGER NPT.NPTS 



IStotus ret by LIO calls 
! LIO dev ice ID 

INumber of bytes of doto reod 
INumber of points to somple 



Declare data buffers 



INTEGER»2 RAW_DATA(5000) 15000 point raw data buffer 

REAL VOLTAGES (5000) 15000 voltages output 



Get o raw_dota buffer of 'npt' values | 

Uses L IOJREAD to reod the A/D 
Note length of buffer is in bytes 

As is the returned data length j 



SYS STAT - LIO$READ(DEV ID. RAW DATA. NPTS. DATA.LENGTH, ) 
I F ( . NOT . SYS.STAT) CALL LIB$SIGNAL(XVAL(SYS_STAT) ) 



Convert the raw dato to voltages using LSP$FORUAT_TRANSLATE_ADC | 



CALL LSP$FORMAT_TRANSLATE_ADC(RAW_DATA, VOLTAGES. NPT. . . ) 



RETURN 

END 
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ooo ooo ooo ooooooooo 



subroutine DETACH(CLK_ID.OEV_ID) 



Subrout in* written by Moj Richard R Ryl**, NPS Student 28 Apr 87 

Updot* 16 j un 87 

This subroutine is intended to detach the internal clock and A— 0 board 
from the main progron and allow them to be used for others application. 

This subroutine should be used with subroutine SETUP 



INCLUDE ’SYSSLIBRARY : LIOSET . FOR' 



Declare local variables 



INTEGER SYS STAT IStatus ret by LIO calls 

INTEGER DEV_ID ,CLK_ID ! LIO device ID 



Detach from the Clock 



SYS.STAT = LIOSDETACH (CLK_ID.) 

I F ( . NOT . SYS.STAT) CALL LIBSSIGNAL(%VAL(SYS_STAT) ) 



Detach from the A/D 



SYS.STAT = LlO$DETACH(DEV_ID. ) 

I F( . NOT . SYS.STAT) CALL LIB$SIGNAL(XVAL(SYS_STAT) ) 
WRITE(«.1) 

1 FORMAT (/I x , 'Both clock and A— D device are detached’/) 

RETURN 
END 
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